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Abstract 

\ This paper elaborates on the achievable schemes of interference alignment in constant channels via finite amplify-and-forward 

(AF) relays. Consider K sources communicating with K destinations without direct links besides the relay connections. The total 
number of relays is finite. The objective is to achieve interference alignment for all user pairs to obtain half of their interference-free 

1 degrees of freedom. In general, two strategies are employed: coding at the edge and coding in the middle, in which relays show 

different roles. The contributions are that two fundamental and critical elements are captured to enable interference alignment in 
^ ■ this network: channel randomness or relativity; subspace dimension suppression. 

I. Introduction 

Since interference alignment (IA) is a new multiplexing gain maximizing technique |fl] and relay is considered a cost- 
T — I \ effective solution for coverage extension and capacity enhancement, recently many researchers have been investigating schemes 
to combine these two techniques 0, 0. Various scenarios have been intensively studied, e.g. decode-and-forward relays 0, 
0, multi-user broadcast channel 0, clustered relays Q, two-way relay selection 0. Among all the scenarios, the model of 
■ multi-user peer-peer two-hop interference channel with amplify-and-forward relays without direct links has received increasing 
^ \ attention due to its wide application in practice (5), 11101 , IPTD . fl~2), which show that relays help the system to obtain higher 
O ■ degrees of freedom (DoF) in high signal-to-noise ratio (SNR) regime. 

While peer-peer multi-user interference channel has been traditionally regarded as a challenging scenario due to channel 
. inseparability iTOl . not to mention with relays in between. Furthermore, current relay networks are categorized into two main 
£>• ' sets: one is that relays are auxiliary links in addition to direct links between both ends 0, lfl4l . fl5l . Ifl6l . 1171 : the other 
(N. is that relays are the only connection path for end nodes without direct links [3 |. In general, relay networks with direct links 
are well structured to construct IA. However, this work considers relay networks without direct links where relay generated 
equivalent channels are quite complicated and poorly structured 0, lfT8l . lfl9l . In that case, many approaches are proposed, 
e.g. (3 | requires an infinite number of relays to eliminate interference; [18| exploits ergodic nature in fading channels; |fT9l 
illustrates interference neutralization in a special 2x2x2 network; [9| uses mean squared error (MSE) numerical method to 
minimize interference. 

Three important features are highlighted in the scenario of this work: the number of relays is finite so that it is impossible 
\ to have the unpractical solution in to eliminate interference with infinite number of relays; each node could only have 
• i-H ■ single antenna to achieve IA; time-invariant channels are also available to achieve IA even with single antenna nodes [1|. This 
" work borrows the two strategies generalized in multiple unicast problem ll20l . which are coding at the edge and coding in 
the middle. In the first strategy of coding at the edge, relays randomly construct equivalent channels while end nodes proceed 
with conventional interference alignment schemes. Compared with other research, the most unfavorable conditions are set in 
this work: all end and relay nodes are single antenna; all channels are time-invariant; relay number is set to be finite, e.g one 
or two. Max-flow-min-cut theorem is not directly applicable in this scenario. In the second strategy of coding in the middle, 
optimization techniques are applied to numerically approach interference alignment with all nodes set to be multi-antenna. 
Our contributions are in two folds: 

A new fact is unveiled that when the network has only one relay, Cadambe-Jafar scheme is not applicable in signal 
vector alignment; Motahari-Khandani scheme is not applicable in signal level alignment; asymmetric complex signaling is 
not applicable in phase alignment. The reason is that relay-emulated channels lose randomness or relativity. Then the problem 
is settled by two ways. On one hand, at least two relays are necessary to emulate qualified channels to generate randomness. 
On the other hand, to generate randomness still in the one-relay channel, space-time type of precoding methods are applicable 
at the edge; the conceptual idea of blind interference alignment is also applicable to fluctuate channels with the only one 
relay. By all these analysis and schemes, a novel unique role of relays in the network is revealed to overcome the channel 
randomness issue. 

A novel solution is proposed for IA via relay coding based on a non-trivial application of existing rank constrained rank 
minimization (RCRM) method lETl . Il22l . Conventional optimization and numerical algorithms are non-convex or hard for 
interference alignment ll23l . 0, 0. Since RCRM method is originally not for relay networks, the convexity is proved for 



(N 



the new application. This novel solution has three advantages: a) it considers DoF at high SNR by subspace dimension 
minimization while other numerical methods could only consider sum rate and mean squared error (MSE); b) it is universal 
to be conveniently applied to a number of scenarios, e.g. MIMO amplify-and-forward(AF) relay channel, MIMO two-way 
AF relay channel, MIMO Y channel, and MIMO multi-hop relay channel, c) actually conventional analytical and numerical 
solutions are hard to obtain for the mentioned scenarios, while this novel method accomplishes the design and is robust to 
poor conditions. Numerical results show its effectiveness. 

II. Problem Statement and Model Description 

A. Basic Model 

Consider K sources and K destinations connected by L relays as shown in Fig. Q] Direct links between end users are 
not available. All nodes are single antenna, denote the channel from source i to relay I as and the channel from relay 
I to destination j as iji. Denote the sets JC — {1,2,..., -fT}, C = {1, 2, . . . , L}, so that i,j £ /C, I € C. All channels are 
quasi-static, i.e. h;; and iji are scalar constants. 
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Fig. 1. Basic model of sources, relays, and destinations 

Time-extended MIMO scheme in [3] is used so that T consecutive symbols form one signal and the relay coding matrix at 
l-th relay is 



Gi = 



i(p, q) 



TxT 



(1) 



in which gi(p, q) is the relay gain from the q-th time slot to the p-the time slot at the l-th relay. Then the received signal at 
j-th user can be written as 



K 



L 

9ji = ^ ijiGihu 



(2) 



(3) 



where X; e C is the d input streams at i-th source and Vj G C is precoding matrix. Zj G C is the noise at 
j-th receiver, and Yj G C Txl . The equivalent channel matrix from the i-th source to the j-th destination is defined as 8^ for 



simplicity. The holistic system equation is also written as: 
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B. General Model 

For more general cases, the model is illustrated in Fig. [2] Denote K nodes on one end as Ii,Is, ■ ■ ■ ,Ik',L relay nodes in the 
middle as Ri, i? 2 , • • • , Rl, K nodes on the other end as J%, ■ ■ ■ , Jr-i, Jk- Define sets JC = {1, 2, . . . , K}, C = {1, 2, . . . , L}. 
All channels are time-invariant. Assume all nodes have all channel knowledge to cooperate. 




Fig. 2. General model of sources, relays, destinations 

In the case when all nodes are single antenna, denote h..R,j 4 , hj-.jj,, as channels from Ii to Ri, Ri to Jj, respectively, where 
i,j 6 JC, I G C, all channel coefficients are scalars. While in the case when all nodes are multiple antenna, denote Hr,^, 
H/^,, Hj^,, Hfl,^, as channels from Ii to Ri, Ri to Ii, Ri to Jj, Jj to Ri, respectively, where i,j 6 /C, I S C, all channel 
coefficients are matrices. 

III. Characteristics of Amplify- and-Forward and Finite Relays 

Before detailed design and analysis, this relay scenario needs to be characterized and some critical features should be exposed. 
First, the applicability of max-flow-min-cut theorem in this network is actually in doubt along with the amplify-and-forward 
scheme. Second, the number of relays, achievable DoF, and coding strategy are correlated in a complicated form. 

A. Amplify -and-Forward Scheme and Max-Flow -Min-Cut Theorem 

A natural question arises here: Does max-flow-min-cut theorem still work in this IA relay network? Another similar question 
is: in a X-pair network with finite relays, what does the achievable DoF equal to? The answers to these two questions are 
quite nuanced actually. 

Look into the K -pair network with L relays. The first hop is a K x L MIMO link and the second hop is a L x K MIMO 
link. Conventionally, both the first and the second channel matrices only have a rank of L, which means that for the whole 
if -pair network, the minimum cut is L. However, the degrees of freedom for the network is not necessarily L. Consider the 
K -pair network regardless of the middle connections, so that there are potentially K DoF for this equivalent network as long 
as with joint processing across the K transmitters or receivers. The critical difference is that the intermediate relay nodes do 
not necessarily process data flows directly, and could act only as equivalent channels. In addition, instead of joint processing, 
a distributed manner such as interference alignment could still possibly achieve K/2 DoF. However, it is also important to 
highlight that when there are finite single-antenna relays, the equivalent channel may be rank deficient, which is discussed 
detailedly in the following parts. 

DoF as in the signal processing model is not simply equal to capacity flow links as in network coding model. In the flow 
network model, each node has separate inputs and outputs, and each edge represents a separate link. While in amplify-and- 
forward relay networks, arbitrary flows could overlap through same relays, and links in each hop does not represent separate 



flows. Max-flow min-cut theorem is originally applied to wired networks with single-letter characterizations. A recent model for 
wireless relay networks is known as the linear deterministic relay network model with a max-flow min-cut result pertaining to it 
11241 . The algorithmic framework is introduced by Avestimehr, Diggavi and Tse, incorporating the key features of broadcasting 
and superposition. The signals are elements of a finite field and the interactions between the signals are assumed to be linear. 
This model is based on linking systems and the max-flow min-cut theorem is applicable with matroid intersection or partition. 

Moreover, [20| shows that in a network with 3 unicast sessions each with min-cut of 1, whenever network alignment can 
achieve rate of 1/2 per session, there exists an alternative approach including routing, packing butterflies, random linear network 
coding, or other network coding strategies instead of alignment. However when there are more than 3 sessions, alignment is 
required to obtain the maximum rate as 1/2 the min-cut, while no other method can achieve it. 

In summary, the max-flow min-cut theorem is not exactly available for this amplify-and-forward system, and DoF is not 
determined directly based on the number of relays. So that it is non-trivial to investigate new schemes in the amplify-and- 
forward strategy. Before further understanding and analyzing the DoF of this network, different situations need to be classified 
as following. 

B. Finite/Infinite Relays, DoF Limits, and Strategies 

For the X-pair amplify-and-forward relay network, all the cases are roughly classified into three categories: when infinite 
relays are provided, full K DoF is achievable, by using relay coding; when finite relays are provided, only fractional DoF is 
achievable as the number of users K grows, by using conventional MIMO precoding approach; when a specific range of finite 
relays are provided, K/2 DoF is possible to be obtained, by using asymptotic interference alignment precoding approach. The 
details are shown in Table Q] and discussed as following. 
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TABLE I 

Classifications according to relay number, DoF Limits, and Strategies 



1 ) Infinite Relays, Full DoF, and Relay Coding: Consider the single antenna if -pair relay network in Fig. Q] with the system 
equation of (0 and the equivalent channel in (01. Let the interference at any receiver to be zero, and then the condition for 
j-th user should satisfy the following: 

Qji = Vi=fj 

L (5) 
i.e. J2^ G l h li = ° Vi^J 
i=i 

Equation (O represents K(K— 1) matrix equations, equalling to T 2 K(K — 1) linear equations. The matrices Gi, G2, • • • , Gl 
contain LT 2 variables. Generally, the equations are solvable when: LT 2 > T 2 K(K — 1), i.e. L > K(K — 1). It could be also 
interpreted as for each set of L elements of {Gi(7i, £2), G2(h, h), • • • , Gl(Ii, h)} there are K(K — 1) equations. Therefore 
in this case, the network achieves K DoF by only using the linear relay coding. 

2) Finite Relays, Fractional DoF, and Precoding: Consider the same network in Fig. Q] However, there're only finite relays 
in this case. Instead of relay coding as above, a scheme of precoding on the equivalent channels is proposed, i.e. following 
the coding at the edge strategy ll20l . Relays generate equivalent channels, and the K transmitters and K receivers only see 
the equivalent channels regardless of relays. Then conventional interference alignment methods could be used directly for all 
K users Q, 11231 . 

The relay gain matrices of (HJ are randomly chosen to be in a generic form (full elements), so that G; and corresponding Qji 
are equivalent to general MIMO channels. The leakage minimization algorithm of [23 1 is applicable to approach interference 
alignment. However, according to the feasibility condition in 11251 for symmetric MIMO channels, the equivalent system in (O 
must satisfy 

T + T- d(K + 1) > (6) 

So that each user obtains DoF bounded by y = -j^i - Numerical results are shown in Fig. [3] in which three cases denoted by 
'Generic' use the leakage minimization algorithm in [23 1 to achieve IA for different numbers of relays: relay number L = 1, 
L = 5 and no relay. Total user number K = 5, and the time extension length T = 21, and each user has d = 5 streams. A 



MIMO network with 21 antennas for all nodes is set as a reference case. When the channels are generated by 1 relay or 5 
relays, each user could obtain the same number of DoF as the reference MIMO channel case roughly as 2/ (K +1). 

3) Specific Finite Relays, Half DoF, and Precoding: Although the feasibility condition in [25 1 limits the DoF in the constant 
MIMO channels, however there is still chance to achieve K/2 DoF for the network actually. The reason is that the relays are 
capable of flexibly constructing desired channel structures such as time-variant channels. 

Manipulate relay gain matrices to be of diagonal structures as following: 



G I =Diag{g I (l > l) > g I (2 > 2),». > g I (T > r)} (7) 

where Diag{ } is diagonal function which place all inputs on diagonal line of a matrix output. Then the equivalent channels 
are constructed in diagonal structure as well: 

L 

Q T1 = ^Diag{f j7g/ (1, l)h Ji ,f i ,g,(2,2)hi< ) • • • ,f 3lgl (T,T)h H } (8) 
i=i 

If we only deal with the relay constructed diagonal channels by using the conventional leakage minimization method as 
before, the DoF results would have no improvement as shown in Fig. [3] It shows three cases denoted by 'Diagonal', which 
indicate the same 2/(K + 1) DoF curves as the equivalent generic MIMO channels where K = 5,T = 21, d = 5. Observe 
that the relay number does not affect the algorithmic IA result either. Theoretically, there has been no conclusions so far on 
algorithmic IA feasibility of diagonal channels ll25l . 




Fig. 3. Coding on the edge with leakage minimization and asymptotic design, with different number of relays 

However, if we deal with the relay constructed diagonal channels by using the analytical asymptotic design (Cadambe-Jafar 
scheme) for frequency selective channels 0], each user could approach 1/2 DoF regardless of the number of K users in the 
network. Take the 3-user network as an example to apply Cadambe-Jafar scheme in this relay network. The system equation 
© needs to be modified a little: let Vi 6 C Tx ( n+1 \ V 2 G C Tx ", V 3 G C Txn , T = 2n+ 1 , where n is a positive integer. 
Thus DoF for the three users are not symmetric here and set as £n+i ' 2n+i ' 2n+i res P ec tively. When n is large enough, each 
user could approach 1/2 DoF. Design the aligned interference subspace at each user as following: 



e 12 v 2 = e 13 v 3 ,e 23 v 3 -< e 21 v 1 ,e 3 2V 2 -< e 3 iVi 



(9) 



Then |1| proposed the following analytical design for each precoding matrix to satisfy (|9): 

a = e 12 e 2 - 1 1 e 23 e 32 1 e 31 er3 1 

Vi = [w wA wA 2 • • • wA™l 

1 o (10) 

V 3 = GaiO^VA wA 2 • • • wA n ] 

V 2 = GaiGga^w wA • • • wA"^ 1 ] 

where w = [1 1 ••• l]t £ c( 2n+1 ' xl . In this way, the three-user network could approach 3/2 DoF eventually. It is 
important to highlight that this scheme could be extended to arbitrary A-user case. Then © and dTob are upgraded to a 
more complicated form accordingly The achievable multiplexing gain of the network becomes ' where 

N = (K — 1) (K — 2) — 1, n is non-negative integer. As a result each user could approach 1/2 DoF. 

Numerical results are shown in Fig. [3] In the 3-user network, let n = 10 so that the time extension length T = 21. There 
are four cases: relay number L = 1, L = 2, L = 5, and no relay. In addition, a case of frequency selective channel with 21 
dimensions for all nodes is introduced as a reference. The four cases using the asymptotic design are denoted by 'Asymptotic'. 
When the channels are generated by 2 relays or 5 relays, each user could obtain the same DoF as the reference frequency 
selective channel case. However, when the channels are generated by only 1 relay, each user obtains zero rate, and zero DoF, 
i.e. the interference alignment scheme fails. 

The reason for the failure of the case of L = 1 is investigated as following. For L = 1 relay, (0) becomes: 

(-)„ (11) 
Notice G; is diagonal as in (0, then the asymptotic design in ( TTOb degenerates to the following form: 

a = fiihi 2 {f2ihiiy 1 f2ihi 3 {f 3 ih 12 y 1 f3ihii{fiih 13 y 1 
A = air 

Vi = [w aw a 2 w ■ ■ ■ a n w] (12) 
V 3 = f2ihn(f2ih 13 y 1 [aw a 2 w ■■■ a n w] 
V 2 = /3i/iii(/3i^i2) _1 [ w aw ■■■ a n-1 w] 

where a is scalar coefficient equal to 1, and 1^ is an identity matrix. Observe in ( fT2b the three precoding matrices Vi, V 2 , 
V 3 have only rank 1 because all are linear dependent to be eliminated, so that the asymptotic solution fails due to the 
lack of channel randomness or relativity in (T). which is the core idea of the IA mechanism. Similarly, in the A'-user case, 
the asymptotic solution also fails when there is only 1 relay. Meanwhile recall the numerical result by leakage minimization 
algorithm in equivalent MIMO scheme, the achieved DoF is low but non-zero. 

Comparing with the case of L = 2, the relay generated diagonal channel in ([8]) becomes: 

e ii = / i iG 1 /i li + / J - 2 G 2 /i2i (13) 

It needs to be proved that Oji are not linear dependent to lose channel randomness with the following lemma. 
Lemma 1: Qji and 0„ m are linear independent almost surely for arbitrary non-identical and (n, m). 

Proof: Assume 
3/3i, p 2 eC.0i,fo? 0, such that ftG^ + /3 2 9„ m = 0. 
Expand and group by Gi and G 2 to: 

(filfjlhu + Mnlhlm)Gi + {Plf j2 h 2i + p 2 fn2h 2m )G 2 = . 

Gi and G 2 are arbitrary generic diagonal matrices, therefore almost surely the coefficients are zeros, i.e. 

(filfjlhu + fcfnihim) = 0- (fiifjshvi + /3 2 f n2 h 2m ) = 0. 
Forming the ratio — 4 1 , so we have / 3l , ,t1 ' = p 2 , 2 ' . 

Since fji,hii,f n i,hx m ,fj 2) h 2 i,f n2 ,h 2m are all random independent scalars, the equality fails almost surely. So that (3i, 
j3 2 do not exist almost surely. Then Qji and Q nm are linearly independent almost surely. ■ 

Notice 'almost surely' cases are recognized as successful interference alignment (TJ. Lemma[T]shows that Oi 2 21 1 ,0 23 9 32 1 , 
and 6310^3 are not scaled identity matrices, therefore their total product, A, is not a scaled identity matrix almost surely 
considering channel relativity of these three terms. As a result the asymptotic design in ( fTOb still works here, rather than 
degenerating to the form ( fT2l . In summary, Lemma Q] reveals that two or more relays can generate channel randomness or 
channel relativity in the equivalent channels with problems in (TTTT) and (fT2l successfully prevented. Then the asymptotic design 
in ( [Tol l can be applied successfully. 



IV. Role of Relay in Precoding Scheme in Constant Channels 

The above 3-user case reveals a critical issue for the constant channel to implement interference alignment. Therefore it 
is important to look into general if -user networks, as well as different class of IA schemes. In this section, for the general 
scenarios of relay networks, the model is shown in Fig. [2] : define nodes Ii, I2, ■ ■ ■ ,Ik as sources and Ji, • • • , Jk-i, Jk 
as destinations, whereas each Ik communicates with each Jk via relays Ri,B,2,- - ,Rl- Relays are half-duplex so that the 
transmission procedure consists of two stages where relays forward data in the second stage. Interference alignment is designed 
with the following strategy: relays construct equivalent channels; source and destination nodes proceed precoding and zero- 
forcing, i.e. coding at the edge as ll20l claimed. The objective is still to approach 1/2 due DoF (exclude duplex factor) for every 
user, which is quite difficult and non-trivial for a relay network under quasi-static channel condition. The work is published 
in ||26| and submitted in |27|. 

A. Time-extended Signal Vector Alignment 

The most common approach to implement IA in this network is the time-extended MIMO scheme as in |3|, where T 
consecutive symbols form one signal vector and IA is designed in the vector space. The above example of 3-user network 
as in equations dT) and <[8J exactly constructs the relay gain matrix G; and equivalent channels 9ji with this time-extended 
MIMO approach. Therefore, it is important to extend the design to arbitrary if -user networks, and look into the same issue 
caused by constant channels with single relay in the network. 

Gi and 0^ could be diagonal or generic. If they are generic, it is equivalent to MIMO channel constrained by feasibility 
conditions |25|, so that each user could only approach 2/(K + 1) DoF. If they are diagonal, the channel is equivalent to a 
frequency selective channel in JTJ, [28] , then the asymptotic design of Cadambe-Jafar scheme could be applied. Notice ll28l 
has the same core design structure as (TJ, so that the design of 11281 is illustrated here for a general if -pair relay network as 
in the following form: 

$ ji = ej 1 1 e ji Qile ia vj,iG/c\{i} (14) 



v 3 = i*23 1 n 



vH n 




(15) 



Vi = e^GiaVg Vie£\{l,3} 

where Vj G c Txd i j s the precoding matrix for arbitrary source ^ 1,3, which contains d{ input streams. A = 
{(j, i) I 3,i e K. \ {1}, j ^ i, (j, i) ^ (2, 3)}. 1 T G C Txl is the all one vector, n* is integer. Notice T = + C*^) 

and N=(K- 1){K - 2) - 1. Vi has a dimension of T x V. ( has a dimension of T x ("*^ JV ), Vi + 1. So that each 

user could approach 1/2 DoF when n* is large for arbitrary K. 

1) The case when the network has single relay, L — 1: However, in the special case when L = 1, this scheme does not 
work. The reason is easily shown in the equations (0 and ( fl4l > as following: 

• 1 T = (h. J;i R 1 Gih.R 1 i 1 )~ 1 (h.j j R 1 Gih Rl i i ) 

■ (^j 1 R 1 G 1 h Rl i 1 )~ 1 (la.j lRl G 1 b. Rl i s ) ■ 1 T (16) 
= h J^R 1 h R 1 1 i 1 h JjRi h RiuKhR 1 h R 1 1 i 1 h JiRi h Rih ■ !t 

Apply ( TTol ) to the precoding solutions of ( fl5l ), then we could observe column subspace of Vi, Vi G K. in ( fT~5b collapses to 
one dimensional space parallel to It- It means the signal subspace of Vi degenerates so that all the transmissions fail to obtain 
DoF. The following remark is summarized. 

Remark 1: For the equivalent channels generated by single relay, interference alignment is not feasible to approach 1/2 DoF 
for every user by using existing asymptotic designs of 11281 , iffl . 

2) The case when the network has at least two relays, L > 2: While for the case of L > 2, asymptotic IA scheme could 
be much likely applied to obtain K/2 DoF of this network for arbitrary K. As a primitive investigation, set L = 2. The key 
is to prove that the column subspace of Vi, Vi G K, in ( fl5l l almost surely maintains its rank. First, we need to look at the 



(Pahj j n 1 iiR 1 i i hj n R 1 hii 1 i 1 hj 1 R 1 hii 1 i m + /3i,hj n R 1 hii 1 i rn 'hj j ii 1 hR 1 i 1 hj 1 ii 1 hR 1 i i ) ■ GxIt 
+ [/3a(h-j j R 1 h.R 1 i i hj Tl R 1 hR 1 i 1 hj 1 R 2 hR 2 i m + hj j R 1 hR 1 i i hj 7l R 2 hR 2 i 1 hj 1 R 1 hR 1 i m + 
hj J fl 2 hB 2 / i hj„ij 1 hij 1 / 1 hj 1 ij 1 hjj 1 j m ) + f3b(hj 7l R 1 hR 1 i m hj j R 1 hR 1 i 1 hj 1 R 2 ]iR 2 i i 
+ hj„_R 1 h_R 1 _f m hj j ._R 2 h_R 2 / 1 hj 1 _B 1 h_R 1 j i + hj„_B 2 h_R 2 / m hj i _R 1 hjj 1 / 1 h.j 1 _R 1 hn 1 j i )] • G1G2I-T 
+ [/3a{iij j R 1 hR 1 i i hj n R 2 hR 2 i 1 hj 1 R 2 hR 2 i m + hj j R 2 hR 2 i i hj n R 1 hR 1 i 1 hj 1 R 2 hR 2 i m + (19) 
h jj r 2 hfl 2 / 4 h j„h 2 hij 2 / 1 h j 1 r 1 hR 1 j m ) + fit (h j n r 1 \ir 1 i m hj. r 2 \ir 2 i 1 hj 1 r 2 \ir 2 i i 
+ hj„fl 2 h_R 2 / m hj j _R 1 h_R 1 j 1 hj 1 _B 2 h_R 2 _r i + hj„_R 2 h_R 2 / m hj ;| ._R 2 h_R 2 _f 1 h.j 1 i{ 1 h_R 1 i i )] • G1G2I-T 
+ (Pah-jj R 2 hR 2 i i hj n R 2 hR 2 i 1 hj 1 R 2 hR 2 i m + j3bhj n R 2 \\R 2 i m hj j R 2 \\R 2 i 1 hj 1 R 2 \\R 2 i i ) ■ G%1t 
= T 



core elements < £ > 23 1< I ) :/ i which constitute Vj. The objective is to prove all the <& 2 3^ji terms are linear independent. Then the 
expanded form is as in equation ( TTTb and proved through the following three lemmas: 

<£-i(jj.. _ (hjafliGih^^ + hj 2 R 2 G 2 h R2ll ) ^ 

23 11 (hj^jGih^j/j + h,/ jfl2 G 2 h fl2 / 1 ) 
(hj J R 1 Gih RlIt + h,7 j / ?2 G 2 h/ ?2 / i )(h,/ lfll Gihi ?1 /3 + hj 1 ^ 2 G 2 h/ i2 / 3 ) 
(b-jaii! + hj 2 fl 2 G2hij 2 j 3 )(hj 1 fl 1 Gihij 1 j i + Ii/^G^i/^t) 

Lemma 2: Gf It, G 2 G 2 1t, GiG 2 1t, G 2 1t are linear independent almost surely. 

Proof: Let G^ 1 G 2 = Diag{Ai, • • • , At}- Diag{-} denotes the function to generate a diagonal matrix with all elements 
in the set. Suppose the lemma is false, then (G^ 1 G 2 ) 3 1t, (GJ~ 1 G 2 ) 2 1t, G^ 1 G 2 1t, It are linear dependent. So that there 
exist non-zero ct\, a 2 , 0.3, satisfying the element-wise equations of a\\\ + a 2 A 2 + a^At + a.4 = 0, W = 1, . . . , T. Since 
Gi, G 2 , Ai, ■ • • , At are random generated independent parameters, and rank-3 polynomials could not have T non-identical 
roots almost surely, so it proves the lemma is true almost surely. ■ 
Lemma 3: ^^^ji^-T and $23 $ nm lT are linear independent almost surely for arbitrary non-identical (n,m) G A, 

where A = {(j, i) | j, i G K \ {1}, j + i, ^ (2, 3)}. 
Proof: 

Suppose 3 non-zero (3 a ,f3b satisfying /3 a ^23 }^ji^-T + 23 ^nm^-T — 0t> where Ot G C Tx1 denotes an all-zero vector. 
Derived from (fl4] i. the following equation is obtained: 

P a ejiQ n iQ lm lT + PbQnmQjlQlilT = T (18) 

By using (01 with (fT8l . further we have equation ( fT~9b . 

In (T% , since Gf It, G 2 G 2 1t, GiG 2 1t, G 2 1t are actually four independent T-dimensional bases, then all their coefficients 
must be zero. We obtain another four linear equations about (3 a , (3b, with coefficients composed of h.j j R 1) 'hii 1 i i ,h.j j R sl ,h.R a i i , 
which are all random independent scalar values of the channels. Therefore the solutions of (3 a , (3b are zero almost surely, which 
contradicts the initial non-zero assumption. So that it proves the lemma. 

■ 

Lemma 4: All vectors in {$ 2 3 1< i 5 jilT, (j, i) G .4} are linear independent almost surely, where A — {(j, i) j,i G /C \ {1}, j / i, (j, i) ^ 
Proof: 

The procedure is similar to Lemma [3] Only notice that the set contains N vectors. The new equation corresponding to (fT8l 
has degree of 2N — 1 over 0j,;. Expand to a new equation corresponding to (fT9l l, similarly there are totally 2N T-dimensional 
bases such as G 2A, 1 T , G 2JV_1 G^1 T etc. Then there are 2N linear equations about N variables j3x, . . . , /3jv, which are forced 
to be zero. So that the independency is proved. ■ 

Finally, in addition, consider the linear independency of the exponentiation terms of $23 ^jIt in the precoders V,, Vi G K 
in (TT~5T > - Because the signal dimension T is large enough to afford the number of bases in space, so that all Vi,Vi G K are 
possibly well-conditioned to implement interference alignment. However, it needs further rigorous proof in our future work. 

In summary, the comparison of the cases L = 1 and L > 2 in signal vector space illustrates that in single antenna constant 
relay channels, single relay is not able to provide channel randomness or relativity which is the key to interference alignment, 
while at least two relays are necessary to provide this feature to approach 1/2 DoF (excluding duplex factor) for every user. 

B. Number-Based Signal Level Alignment 

Besides signal vector alignment for single antenna constant channels, there is another major class of schemes called signal 
level interference alignment. There are several approaches to investigate signal level alignment, including lattice coding, 
deterministic models, and number theory. Among these approaches, IA on the number domain is a very novel and canonical 



method. While our focus in the following work is to show, in relay connected constant channels, in a similar manner to vector 
alignment scheme, how level alignment would also face the applicability issue of existing designs when there is only one relay 
and the issue solved when there are at least two relays. 

1 ) Basic Concept and Design Procedure: At first, signal level could be viewed on the rational number scale which represents 
infinite fractional DoF 11291 . Then Khintchine-Groshev theorem also reveals that the field of real numbers is rich enough to 
be equivalent to vector space to design IA. Furthermore, [30 Theorem 7] uses a generalized version of Khintchin-Groshev 
theorem to extend the designs to complex channels. In real channels, DoF is defined as d^ ca \ = limp_ s . 00 i/ 2 ( \og n (P) wn il e in 
complex channel DoF is defined as ^complex = limp->oo T^jtfp) wnere Csum is sum capacity and P is transmission power. 

By using this Motahari-Khandani scheme in (29), 11301 , K-user single antenna constant channels could approach K/2 DoF. 
Since the objective in this work is to study the role of relays, the core structure and procedure of IA are briefly described in 
a complex channel setting. 

Define u^, = [u^u^ ■ ■ ■ ui ] as the number of d k datastreams sent from node Ik to node Jk, k € 1C. G (— Q, Q)z, 1 < 
m < dk, i.e. belongs to an integer constellation. Each datastream is multiplied by a number z^'™' 6 C, which is called a 
modulation pseudo-vector serving as distinct directions. In order to satisfy the power constraint and control the minimum 
distance of the received constellation, transmission signals should be scaled with a constant A. 

In the meanwhile, each relay node R\ generates a random gain to be a rational number g\ G Q and the equivalent channel 
from ij to Jfc is as following: 



L 

hi 



h JM ■ 9i ■ W< ( 2 °) 

l=i 

Then the received signal at destination node Jk is: 

y k = X9 kk J^4 m ' ) u^+ A^X>| ro) <4 m) (2D 

m—l i=l J i^k m=l 



According to 11291 , 11301 , the structure design of pseudo-vectors vY™' on number domain is similar to the design of V, in 

e as well. In this design, all v?"^ 

< Uji < n — 1 i = k,j =/= k 



vector space of (28], however there exists notable difference as well. In this design, all f k , 1 < m < d k belong to a set B 



^=m n 

j = l 



K K 



< ctji < n Otherwise 



(22) 




n is integer, so that the number of streams d k = \B„ k \ — n K ^ 1 {n + l)^^ 1 ) 2 . 

For destination node the received signal space of (fJTJ contains desired signal subspace formed by 9kkB Uk and interference 
subspace formed by 9kiB Vi ,i ^ k. Observe all the interference subspaces overlap to the same set of B' Vk . 

< aji <n\ (23) 

So that \B' Vk \ = (n + l) K ( K ~ x \ Meanwhile the desired signal subspace 9kkB Vk and interference subspace B' Vk are distinct. 
According to 129, Theorem 6], the total DoF is „, Jr " „ — r?rg r; — . which approaches K/2 when n is large 



1 («+l)< K - 1 )^+(n+l) A '< A '- 

2) Issue in the Three-User Special Case: In the three-user case, there is a special form of design in ll29l Definition 1] 
presented as following: 

yi — #i • xx + x 2 + x 3 + zi 
y 2 = &2 ■ %2 + Xl + x 3 + z 2 

2/3 = $3 ■ %3 + xi +tf ■ x 2 + z 3 ( 24 ) 
a 9\ 3 9 2 \9 32 9\\9\ 2 9 23 9 22 9\ 3 9 33 9\ 2 9 2 \ 

= ft ft ft ' 1 = ft ft o = ft ft ' 3 = ft ft ft 

<712<723t'31 fl2f2lfi3 V\ 2 V 23 ^YlV 23 {) 3 \ 

So that user 1 has the transmit directions {1, $o, i?q, . . . , $q }, user 2 has the directions of {1, §q, . . . , $q -1 }, and user 3 has 
the directions of {1, z9 , • ■ • , ^q 1 }- 



(hj 1 R 1 hR 1 i 3 hj 2 R 1 hn 1 i 1 hj 3 R 1 hR 1 i 2 — hj 1 R 1 liR 1 i 2 hj 2 R 1 hR 1 i 3 hj 3 R 1 hR 1 i 1 ) ■ g 1 
+ [(h-.j 1 R 1 hR 1 i 3 hj 2 R 1 hR 1 i 1 hj 3 R 2 hR 2 i 2 + hj 1 R 1 hR 1 i 3 hj 2 R 2 hR 2 i 1 hj 3 R 1 hR 1 i 2 + 
b.j 1 R 2 hR 2 x 3 b.j 2 R 1 b.R 1 x 1 hj 3 R 1 b.R 1 i 2 ) — (hj 1 R 1 h,R 1 i 2 hj 2 R 1 hR 1 x 3 hj 3 R 2 h.R 2 x 1 
+ h.j 1 R 1 \iR 1 i 2 h.j 2 R 2 h.R 2 i 3 \i.j 3 R 1 \iR lIl + hj 1 R 2 hR 2 i 2 hj 2 R 1 hR 1 i 3 hj 3 R 1 h Rl i 1 )] ■ gfg 2 
+ [(hj 1 R 1 hR 1 i 3 hj 2 R 2 hR 2 i 1 hj 3 R 2 hR 2 i 2 + hj 1 R 2 hR 2 z 3 h.j 2 R 1 hR 1 z 1 h.j 3 R 2 h.R 2 i 2 + (26) 
hj 1 R 2 hR 2 i 3 hj 2 R 2 hR 2 i 1 hj 3 R 1 hR 1 i 2 ) - (hj 1 R 1 hR 1 i 2 hj 2 R 2 hR 2 i 3 hj 3 R 2 h R2 i 1 
+ hj 1 R 2 hR 2 i 2 hj 2 R 1 h Rl i 3 hj 3 R 2 hR 2 i 1 + hj 1 R 2 hR 2 i 2 hj 2 R 2 h R2 i 3 hj 3 R 1 hR 1 i 1 )] ■ gigl 
+ (hj 1 R 2 hR 2 x 3 h.j 2 R 2 hR 2 x 1 h.j 3 R 2 hR 2 x 2 - hj 1 R 2 hR 2 i 2 hj 2 R 2 hR 2 i 3 hj 3 R 2 hR 2 i 1 ) ■ gl 
= 



If the network has only one relay, L = 1, according to ( f20b and ( 124-b : 



h,7i gi hj?! j 3 hjggi hj?i /1 hjgBi hfii I 2 ' ff 1 
h./ 1 _R 1 ri fll 7 2 h,7 2fll h fll 73hj 3fll h fll 7 1 • gf 

Then the encoding set {l,$o,$o; • • ■ 7^0} use d by a ^ three users collapses. The scheme of ||29l , |30| is not applicable to 
this case to support normal transmission with interference alignment. 

•If the network has two relays, L = 2, IA scheme is applicable as proved by the following lemma: 

Lemma 5: do ^ 1 holds almost surely, for random independent parameters Yij-r x , hj j R 2 , h^j^ b-R 3 j i; g\, g2, where 
i,j e K. 
Proof: 

Suppose do = 1, from (|20| > and d241 l. we have ^13^21^32 — #i2#23#3i = 0, and expand it to the equation (1261 , 
In equation (1261 . since {gf, g 2 g2, gig 2 j 9%} are independent number bases with measure of one, the corresponding coefficients 
must be all zeros. Observe gf and g\ in ( f26b , their coefficients are indeed 0; while observe g\g% and g\g\, their coefficients 
are non-zero with measure one, because of random selected parameters of hj 7^, hj r 2 , h^j^ 117227;. As a result, it proves 
that i? ^l, 

■ 

In summary, the comparison of the cases L = 1 and L = 2 shows that single relay is not able to provide channel randomness 
or relativity in the signal level alignment on number domain, while two relays are capable of provide this feature to approach 
1/2 due DoF for every user. 

3) Issue in the K-User General Case: Extend the study of the three-user case to the general i-T-user case. 

•If the network has only one relay, L = 1, according to (l20l i and (|22t . the set £>„ fc is presented in an equal form of set B„H 



. 1 J-L ~w~y < an <n — li = k,j^k 



3 = 1 



< oiji < n Otherwise 



(27) 



Counting the cardinality \B^}\ is complicated so that it is natural to think of attain an upperboundby relaxing and expanding 



the set 23" to a new set B%}. B^} C Bfl} 



< 7 < nK 2 , < tj, ft < uif 



c = ( n hT uj • (37) • ( n h w 1 e = 7, e « = ^ (28) 



Notice that B" C B^, so that [S^l < \B^\. Further loose the set by: 

C 7 = {( T i) T 2, • • • ,tk)\0 < Tj < nK,Vj e IC, E T i = 7} 

je/C 

v i = {( t i,t- 2 , . . - , Cjc)|0 < ft < nK,Vi £ JC, q = 7} 

ie/c 

raif 2 nK 2 

< E 1^1 • 1^7 1 < E( n * + l ) K • ( nK + !) J 

7=0 7=0 

Compare with the original due cardinality of > 



(29) 



IB" I K 2 + l)(nK + l) 2g ^ooj^oc 



L (n + 1)(^-!) 2 n^- 1 (n + l)(^" 1 ) 2 



It reveals that in equivalent channels generated by single relay, L = 1, the encoding set B Uk degenerates to negligible size 
for large n and K, i.e. Motahari-Khandani scheme ||29l , ll30l is not applicable in this case to support normal transmission with 
interference alignment. 

•If the network has at least two relays, L > 2, we have similar conclusions as lemma [5] such that the precoding set are 
greatly expanded with independent elements to support normal transmission with interference alignment. It also needs further 
comprehensive validation and rigorous proof with exact requisite conditions. 

C. Asymmetric Complex Signaling Alignment 

Asymmetric complex signaling alignment is categorized into a class of phase alignment. It is proposed by Jafar as an 
effective measure to solve the problem of obtaining high DoF in constant interference channels ll3~Tl . However, the work in 
this thesis investigates if it is applicable to the relay generated equivalent interference channel. 

It has been conjectured by H0st, Madsen and Nosratinia that complex Gaussian interference channels with single antenna 
nodes and constant channel coefficients have only one degree-of-freedom regardless of the number of users [32|. Intuitively, this 
conjecture indicates the optimality of orthogonal medium access such as TDMA where each user is assigned a fraction of DoF 
and the sum DoF is equal to one. Original interference alignment (TJ, ll33l provided a powerful tool to achieve K/2 DoF for 
K-usei interference channels, only under the condition of time-varying or frequency-selective channel coefficients. However, 
in single antenna constant channels, it was not known if H0st-Madsen-Nosratinia conjecture is right until the appearance of 
asymmetric complex signaling method. 

By using asymmetric complex signal inputs, Jafar shows that at least 1.2 DoF are achievable on the complex 3-user 
interference channel with constant channel coefficients PP . without any assumption of time-variations/frequency-selectivity 
used in prior work. In conventional wireless systems only circularly symmetric complex Gaussian random variables are typically 
used in order to maximize entropy, while in asymmetric complex signaling the inputs are chosen to be complex but not circularly 
symmetric. For example, in an specific case of phase alignment, each transmitter uses only real valued Gaussian inputs, and 
ensures that interference at each receiver aligns in the imaginary dimension when the desired signal is in the real dimension 
of the complex space. This idea could be extended to general values of channel coefficients thereafter. 

Detailed introduction is omitted here scheme, while the basic model and conclusion is briefly presented. First, each received 
complex signal is expressed in real values as following (noise is ignored for simplicity): 

Re{Yj} 
Im{YO 

Defined as U{4>ij ) 

where the complex output Yj and complex input Xj are split into scalar counterparts of real values, and the complex scalar 
channel H,j is expressed by its amplitude \R%j\ and a transformation matrix U(cf>), which is a rotation matrix defined with the 
following properties: 

= U{-<t>) 

U(4>)U(9) = U{6)U{<t>) (32) 
U(4>)U(6) = U{4> + 0) 

With the new formulation, it is proved that the 3 user complex Gaussian interference channel with constant coefficients 
achieves 1.2 degrees-of-freedom if all of the following conditions are satisfied ||3T1 Theorem 2]: 

<t>Jxh + 4>Jvh - 4>Jii 3 - <t>Jxh / 
<f>JxH + <t>J 3 h ~ fijsh ~ 4>Jxh j 1 
0./ 2 /i + <t>J 3 h - <Pj s 1i ~ / 

(f>j 2 h + <Pj-lI2 - <f>Jih - 

(f>J 3 I 2 + &J-LI3 - - &J3I3 / 

0./3J1 + <t>,hi 3 - <pj 2 ii - <pj 3 i 3 ^ 

where I\, I2, I3 are transmitter nodes and J±, J2, J3 are receiver nodes. Notice the coditions of d33l are naturally satisfied 
almost surely for general 3-user networks. However, consider the relay networks as in Fig. [2j if there is only one relay, then 
the equivalent channel phase from transmitter Ii to receiver Jj is equal to the sum of phases from transmitter li to relay R\ 
and from relay R\ to receiver Jy. 
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(31) 



mod (tt 
mod (tt 
mod (tt 
mod (tt 
mod (tt 
mod (tt 



(33) 



fyjjli = ^Rxh + fijjRi mod ( n ) 



(34) 



Therefore, check the conditions of (|33l again with d34l ), and it is easy to find out the following new results which make the 
design no longer available: 



(t>Jih + <f>Jih 
fijth + 4>J 3 h 

<t>J 2 Il + <t>J 3 l2 
<t>J 2 l3 + 4>Jll2 

4>.j 3 i 2 + 4>Jii 3 

Similar to the time-extended signal vector scheme and the number-based signal level scheme, the remark for complex phase 
alignment is concluded: 

Remark 2: For single antenna constant channels, if the X-pair users are connected by only one relay, then asymmetric 
complex signaling scheme would not work to obtain more than one DoF in this relay generated equivalent channel. 

V. SOLVEMENT OF SINGLE RELAY ISSUE IN PRECODING SCHEME IN CONSTANT CHANNELS 

Previous analysis shows that single relay is not capable of providing channel randomness to support interference alignment 
schemes by precoding at transmitters. However, in this section, we look into possible alternative solutions in the scenarios 
only with single relay. In order to address this issue, a natural strategy is to transform the role of multiple (L > 2) relays to 
the role of multiple (L > 2) time slots scenario with one relay. 

As a primitive investigation, at first the following part only looks into the case when there are only two separate network 
statuses. In particular, the case is equivalent to either having two relays, i.e. L — 2, or dividing the whole transmission period 
into L = 2 parts. Suppose the transmission period starts from time slot 1 to time slot T (set T as even). To divide all the 
T time slots into L = 2 parts, just think of the channel statuses at odd time slots as one virtual network with one virtual 
relay; then think of the channel statuses at even time slots as another virtual network with another virtual relay; add these 
two classes of statuses correspondingly to form a virtual network with two relays. Then the equivalent channel as shown in 
previous equation (0 is presented as: 

%= Diag{fWg l (l^)h|Jl,fj?]g l (3,3)h|f,... ,f$- 1] gl (T-l,T-l)h%- 1] } 
+Diag{f^gK2,2)h^^ 

In d36*l i, actually I = 1 because there is actually one relay, hj*' represents the channel from source node i to relay node I at 
the time slot t, and f|| represents the channel from relay node I to destination node j at the time slot t. 

[11 f2l \T] 

Notice that an issue arises here. Because the channels are constant in all T time slots, so that = = • ■ ■ = and 
ft, = v3 = ■ ■ ■ = ft, . Then the above equivalent channel ( f36t becomes: 

% = fjJ ] h{J ] Diag{g,(l, l),gi(3, 3), • • • >gI (T - 1,T - 1)} 

+fjM 1 Diag{g I (2,2),g l (4,4),... ,gi(T,T)} (37) 
^ f|.; 1 hL 11 Diag{ 3 (l), ff (2),..., 3 (r/2)} 

Observe the final combined result in ( |37| i and find that the channel is equivalent to previous case of ( fTTT i which has only 
single relay. In this case, it has been proved interference alignment is definitely not applicable at all. 

Therefore, to artificially emulate the case of two relays where the channels should show randomness, two solutions are 
proposed in the following two parts. One is by using switchable antenna for the relay node in the middle, so that each time 
after the antenna is switched all the corresponding channels change together in a random manner; the other is to proceed 
fluctuation precoding scheme for both end nodes at the edge with double-layered symbol extensions, so that each time the 
precoding gains change all the channels change randomly as well. 

A. Antenna-Switching Coding in the Middle 

The solution of antenna-switching coding is inspired by the idea of reconfigurable antenna and blind interference alignment. 
Reconfigurable antenna is an emerging technology, changing its characteristics by switching geometrical-metallic segments 
to create pre-determined independent modes in every time slot |34l . Each distinct configuration is a different mode. This 
technology is recently introduced to the research of blind interference alignment (35), which utilizes reconfigurable antenna to 
manipulate the channel directly to create channel fluctuation patterns that are exploited by the transmitter. 
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(39) 



Therefore, the relay Ri (I — 1) is set as a reconfigurable antenna and it has two modes which operate for odd time slots 
and even time slots respectively. An important feature is that channels change between odd time slots and even time slots, 
i.e. h™ f hg> and f + if while = if =... = faM if = if = = if and if = ff = . • • = and 
fjj' = fj 4 ' = •• • = fjp. Then the above equivalent channel d36l l becomes: 

©ii= fj} ] h{J lDia g{s(l»l).gi(3,3),---,gi(T-l,T-l)} 

+fjf ] rf Diag{ gi (2, 2), g,(4, 4), ■•• , g,(T,T)} 

Since it could not be further combined, so that it is equivalent to the case with two relays in equation ( TT3l . Interference 
alignment is much likely to be applied in this case. However, regarding DoF, there is an additional factor 1/2 due to the 
double-layered symbol extension. Thus the X-pair network could approach K/A (excluding duplex factor) DoF with one 
single antenna relay on condition that IA is finally successful. Notice the DoF is still scalable regarding K. 

B. Antenna-Fluctuating Coding at the Edge 

The above antenna-switching solution has an extra requirement of the antenna being physically reconfigured. While in this 
part, a simple but effective solution is proposed to directly proceed coding at the end nodes, instead of the relay. Let each 
source node i has an additional random gain of of 1 at the time slot t; and each destination node j has an additional random 
gain of /jj 4 ' at the time slot t. Then the equivalent channel (l36l > becomes: 

% = 1^^(1,1)^^^(3,3)^, 

-1,T -!)}$-*<$■-*} 

+Diag{/f]ffgK2,2)hl?af^]fWgK4 ) 4)hWaM ) 

-,^a(T,T)W} 

Similarly, let af ? af ] and fif * fif, while af = af =■■■ af af = af = • • • af 1 and = ff =■■■ = 
P? ] = /S' 41 = •• • = Besides, because channels are constant, h\] ] = hj? 1 = • ■ • = hjf 1 and = if} = ■■■ = if} ] . 

J 2 J J 1% 1% 1% JL JL JL 

So that d39l becomes: 

%= ^ 1 ]fWhWa| 1] Diag{g I (l,l) s g I (3,3),...,g I (r-l,T-l)} 

(40) 

+Pf ] ifh§af ] Diag{ a (2, 2), gI (4,4), • • • , gi(T, T)} 

Then it is also equivalent to the case with two relays in equation ( f]~3b . Interference alignment is much likely to be applied 
in this case. Also, the JT-pair network could only approach K/A (excluding duplex factor) DoF which is scalable regarding 
K on condition IA is finally successful. Notice that in this solution reconfigurable antenna is not required. Simple coding at 
the edge could generate channel randomness to apply IA schemes. 

VI. Numerical Relay Coding Scheme aiming at DoF 

Although the analytical design could approach K/2 DoF, it is very complicated to implement and requires very large 
dimensions in signal subspaces so that it is not quite practical. Thus, in the following part, a numerical method is proposed 
for the relay coding, corresponding to the strategy of coding in the middle ll20l . 

Relay coding and optimization methods have been studied in a lot of works. However, all the methods take rates as direct 
objectives. While concerning the high SNR region, in term of DoF, a novel idea is naturally considered to optimize the achievable 
DoF directly with relay coding. However, it is necessary to find an appropriate tool to support the design. Therefore, in this 
work, it is proposed to utilize a powerful tool called rank constrained rank minimization (RCRM) ETTl . 11221 . While RCRM 
method is originally used to design the precoding scheme to achieve interference alignment in general K-usei interference 
channels. So the previous solution is not applicable to the relay coding problem in our scenario. 

In this work, the contributions are two-fold: first, the RCRM formulation ||2~T1| is transformed to a novel problem d47b 
regarding new decision variables of G;, and it is proved that its convexity is maintained to proceed numerical algorithms; 
second, along with a variety of extensions of this scenario, a general relay coding algorithm is presented to solve such kind 
of problems in a universal practicable way. 



A. Interference Alignment by Rank Constrained Rank Minimization (RCRM) 

Consider the model in Fig. Q] with all nodes equipped with single antenna. K sources communicate with K destinations via 
L relays. Time-extended MIMO scheme is used with a period time of T. The relay coding matrix is G; as in ([T). Then the 
equivalent channel Qji is shown in (01, composed of G;. Each user has d streams in the scheme, so that Vj G C Txd is the 
precoding matrix at source i; Gi G C TxT ; Uj G C Txd is the receive filter matrix at destination j. Then each receiver j linearly 
processes the received signal by zero-forcing. The received interference and signal should satisfy the following conditions: 

+ (41) 

rank(Ute i3 -Vj) = d 

Define interference and signal subspace respectively for receiver j: 

where Jj({Vi}|Sj 5 {Gi}^ =1 ,Uj) gathers all interference observed at the j-th destination and Sj(Vj, {Gl}f—i, Uj) gathers 
the desired signal from the j-th source. However, compared with original problem in Ell , notice in our work, the formulation 
(l42l i contains G; as new variables through the equivalent channel Qji as in (0, in addition to prior variables Vj and Uj. 

The conventional method of IA to approach zero-forcing conditions in ( HTb is the interference leakage minimization algorithm 
l23l . While the idea of |211 . 1221 indicates that, instead of leakage minimization, the conditions of (|4"TT > could be dealt with 
directly as a rank problem as following: 

rank(^({V,}^., {G;}f =1 ,Uj)) = 
rank(5 7 (V j ,{G i }f =1 ,U,)) = rf 
Then the problem ( |43T > is further formulated as a rank constrained rank minimization as following: 

K 



minimize rank( t 7 ) 



(44) 

s.t. : rank( t Sj) = d, Vj G K 

The problem (l44l guarantees that the interference subspaces collapse to the smallest possible dimension, under the constraints 
of the desired signal subspaces to have full ranks. While minimizing the sum of ranks of the interference matrices is equivalent 
to maximizing the sum desired DoF of the network. 

Since the objection function in (144-b is nonconvex, ||2TI proposes a tight convex approximation regarding J.y. 



Erav(jSank(.3)) =^||^||* 

^ K d 

-EE"-**) 



f 1 3 = 1 n=l 

where conv(-) denotes convex envelope of a function, and || ■ ||* = 2„ = icrn(') is tne nuclear form of a matrix, which 
accounts for the sum of the largest d singular values, /i is a constant satisfying [i > max je ^ &i{Jj)- 

Since the constraints in d44l) is also non-convex feasible set, lETl provides an approximation of convex set regarding Sj : 

Amin(iSj) > £ 

9 >- n (46) 

where A m i n (5j) is the minimum eigenvalue of Sj and e is a small positive constant, and Sj >z Odxd denotes the matrix is 
hermitian positive semidefinite. 

With the above convex approximations (l45l l and d46l >. the non-convex problem in (l44t is relaxed to a total convex optimization 
ll36l . d45b suppresses interference spaces to the smallest dimension, while (|46| | guarantees desired signal subspaces. It is 



equivalent to maximizing DoF of the network. Then the original coding design problem is formulated as a rank constrained 
rank minimization problem as following: 

K 

minimize >J \\Jj\\* 

i =1 (47) 
s.t. : A m i n (iSj) > e, 

Sj h Odxd, Vj e K. 

Notice that the optimization problem of ( |47] i is convex regarding the variables Jj and Sj . In the general IA problems as in 
II2TI . it is also convex regarding variables Vj and U,-. While for relay networks, it is not still guaranteed if it is convex for the 
relay coding variables Qji and G/. 

B. Design of Novel Relay Coding with RCRM in Time-Extended Single Antenna Channels 

Compared with [21], hereby we specify three groups of decision variables G;, Vj, in the optimization problem d47l i. In 
the following part, ( l47| > is to be proved to be convex regarding decision variables G;, as well as original convexity of Vj and 
Uj already indicated in |2T|. The key of the proof relies on the affine relationship between new variables and the objective 
function and constraints, which are shown in lemma [6] and lemma Q respectively as following. 

Lemma 6: The objective function d43T > is convex regarding variables {G/}^_ 1 . 

Proof: 11361 §3.2.2] indicates the operation of composition with an affine mapping preserves convexity or concavity of 
functions. Suppose a function / : M™ — > R, and a matrix A £ M. nxm , and a vector b G R". Define a function g : R m — > R by 
g(x) = f(Ax + b), with domj = {xjAa; + b G dom/}, where x G R™ and dom denotes the domain. Then if / is convex, 
so is g. 

Extract all T 2 L elements from all G; as Q G C T Lxl , and extract all d(K — l)dK elements from all Jj as one vector 
£ G C d2( - K ~^ Kxl . According to © and (g2]i, we obtain that £ = AG, where A G £ d2 ^ K - 1 ) KxT ' 2L is composed of elements 
from {Uj}, {Vi}, {fji}, {hu}. Since £, A and Q are in the complex domain, they could be transformed into £g, As. 
and Gr in the real domain by simple matrix operations. Then £r = ArGr, where g C 2T ixl , £r G C 2d (k- 1 ) 1 ^* 1 ^ 
Ar G c 2d ( K -V Kx2T L . Thus, £m is an affine mapping from Qtr- ETI concludes d45l l is convex regarding £r, therefore (|45] > 
is also convex regarding Qr, i.e. regarding {Gi}^_ x . ■ 

Lemma 7: The feasible set ( |46b is convex regarding variables {Gi}f =1 . 

Proof: ll36l §2.3] indicates that affine function preserves convexity of sets, or allows us to construct convex sets from 
others. S C R" is a convex set. If the function / : R k -> R n is affine, i.e. f(x) = Ax + b, where A G R nxk , x G M fcxl and 
b G R™, then the inverse image of S under /: = {x\f(x) G 5}, is convex. 

Similar to lemma® extract elements from {Sj}, {G;} as 5r, CJr in the real domain. According to (O and (l42l . 5r = Sr^r, 
where Br is composed of elements from {Uj}, {V^}, {fji}, {hu}. II2TI concludes that feasible set of Sr under constraint 
(|46*T i is convex, so that Qr feasible set as the inverse image of Sr is also convex, i.e. {Gi}f =1 feasible set is convex. ■ 

1) Numerical Results for the Single Antenna Network: Thus lemma [8] and lemma [9] prove that (l47l i is a typical convex 
optimization regarding {Gi}[ =1 . So that it is convenient to use the software tool CVX to proceed the numerical results ll37ll . To 
compare the effects of precoding at the edge and relay coding in the middle, 8 different cases are simulated with configurations 
shown in Table HI] 
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TABLE II 

Configurations of different simulation cases for relay coding 



In Table [HI five different processes are provided in the comparison as well. Process 1 is set as a reference, where the relay 
coding matrices {G;}, precoding and zero-forcing matrices {Vi} and {Uj} are randomly generated. While in process 2, only 
{Vi} and {Uj} are randomly generated, relay coding matrices G; are designed by solving the RCRM problem d4Tb with CVX 
ll37l . Process 3 is on the contrary, where {G;} are randomly generated, while V,; and Uj are solved by using the conventional 
leakage minimization algorithm of interference alignment [23 1, on the equivalent channels Qji as in (0. Process 4 performs 
the leakage minimization algorithm to optimize V, and U 3 additionally on the basis of Process 2. Process 5 applies the RCRM 
method to optimize relay coding {G;} additionally on the basis of Process 3. 
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Fig. 4. Numerical Results of the RCRM method in Time-Extended Single Antenna Channels (8 cases configured as in Table ITU 

Numerical results for the above 8 cases in Table HTl are illustrated in Fig.|U Case 1 shows an equivalent interference channel 
as in (fJJ without any interference alignment design, so that the sum rate deteriorates in high SNR and the achieved DoF is 0. 
Case 2 shows that process 2 uses the rank minimization method to design {G/} to achieve IA and obtain DoF. Case 3 shows 
that process 3 uses the leakage minimization method to design Vi, Uj to achieves IA as well, but obtaining lower DoF. Case 4 
attempts to additionally optimize Vi, Uj by leakage minimization on the basis of process 2, however the DoF is reduced to the 
same level of process 3. The reason is the IA structure formed by relay coding using rank minimization is destroyed, due to 
the inseparable inter-relationships of relay coding variables as shown in (|42T> . Case 5 shows that process 5 takes an additional 
step to optimize {G/} using rank minimization method on the basis of process 3. Notice the DoF rises to the same level of 
process 2, which means the step of relay coding plays the key role in achieving interference alignment. Finally, compare the 
cases 6, 7 and 8 with case 2. Case 6 shows that when the number of users increase from 3 to 4, interference alignment is no 
longer feasible. So that less concurrent streams are allowed, i.e. d is reduced from 4 to d = 2 to resume IA. Case 7 shows that 
when the number of users increase to 4, while maintaining the same number of streams as case 2, IA could be still feasible if 
we increase the time extension length from T = 10 to T = 15. However the DoF is yet not increased along with the number 
of streams, because time extension should be considered resulting in DoF = d/T. So that in case 8, in order to obtain the 
same number of DoF per user as in case 2, the issue is settled by increasing the number of relays from L = 5 to L = 10. 

C. Applying Novel Relay Coding to Diverse General MIMO Relay Channels 

The novel application of rank constrained rank minimization (RCRM) method is proved to be very effective in the relay 
coding of the A'-pair time-extended single antenna networks. Furthermore, it is natural to think of extending this idea to 
more general and realistic scenarios. Network with all nodes equipped with multiple antennas (i.e. MIMO networks) are 
usually investigated in other works |9j, Q. It is in general hard to design beamforming to mitigate and align interference 
with conventional optimization techniques. For example, the end-to-end sum-rate maximization problem is non-convex and 
NP-hard, and implicit MSE(minimum square error) metric and the interference pricing method are considered to obtain sub- 
optimal solutions ID, Q. RCRM method is a promising efficient approach to deal with such MIMO interference channels 
Ell . l22l . It is yet non-trivial to apply the RCRM method to MIMO interference channels with relays, where the convexity for 
relay coding is an issue. Anyhow, the novel approach to apply RCRM in MIMO relay networks provides an simple accessible 
way. 

The following work focuses on general MIMO relay networks without any time symbol extensions. Although it is usually 
infeasible for MIMO networks to obtain 1/2 DoF for every user by interference alignment [25], the idea of suppressing 
interference subspace dimension could be still used to obtain high DoF, which captures the key value of interference alignment 
at high SNR. While RCRM method provides such a measure to minimize the subspace dimension to obtain suboptimal DoF 
through convex approximation 1211 . |22l . The RCRM method is non-trivially extended in this work to apply to relay-connected 



MIMO networks. So that relay coding is involved and designed along with precoding at source and destination nodes, which 
reflects the strategy of coding in the middle l20l . 

Moreover, this novel design is further found to be naturally applicable to a wide range of typical relay networks, which are 
originally hard to analyze and solve. These scenarios include MIMO relay channel, MIMO two-way relay channel, MIMO Y 
channel, and MIMO multi-hop relay channel. 

1 ) MIMO Relay Channel: MIMO relay channel exactly means a two-hop interference channel with half-duplex relays and 
a two-stage transmission procedure. Consider the model in Fig. [2] There are K sources denoted from I\ to Ik', K destinations 
denoted from J\ to Jk', L relays denoted from R\ to Rl- All nodes are set to be equipped with multiple antennas. has 
antennas, Vi G JC; Ri has M x antennas, VZ G C; Jj has M w antennas, Vj G K.. All the channels are constant. H Ri l an d Hx,fl ( 
as channels from Jj to Ri and from Ri to Jj respectively. Hr,^ G C M * xM? and Hj^r, G C m "> xM *. Then the equivalent 
channel 0« G C M » xM i is expressed as: 

L 

e ii = 5^Hj jJll G,H Jl|/l (48) 

2 = 1 

For each user k, Ik transmits d streams to Jk- So that set Vj G <C A k xd as the precoding matrix at source If, G; G C M x xM * 
as the network coding matrix at relay Rf IL, G C M " xd as the receive filter matrix at destination Jj. 

The procedure to apply RCRM method in this MIMO relay network is the same as above in the time-extended single antenna 
network. The objective function and feasible set are in the same expressions of (l45l l and d46b . The final problem is also in the 
expression of d47l >. It is necessary to prove the convexity of the problem in MIMO relay networks by following lemma [8] and 
lemma [9] 

Lemma 8: The objective function d45l > for MIMO relay networks is convex regarding {G;}/^. 

Proof: |[36l §3.2.2] indicates the operation of composition with an affine mapping preserves convexity or concavity of 
functions. Suppose / : M" -> R, A G R ,ixm , and b G R n . Define g : R m -> 1 as g(x) = f(Ax + b), with domg = 
{x|Ae + b G dom/}. Then if / is convex, so is g. 

Extract all M X 2 L elements from all G; as Q G C M * 2 Lxl ; extract all d(K-l)dK elements from all Jj as T G ^(k-^kxi ^ 
According to (gHJ and (02), obtain that T = WG, where W G £ d2 ( K ^) K * M x 2 L [ s composed of elements from {U,}, {V 4 }, 
{Hj flj}, {Hfl,/;}. Express T, W, Q in complex domain to equivalent JF^, Wr (transformation via matrix operations), in 
real domain, so that Jr = WrOr, where G M G C 2M ^ Lxl , Jr G C^ 2 (k-i)Kki ( Wr £ C 2d 2 (K-i)Kx2M x 2 L_ So j- r is an 
affine mapping from C?r. ||2T1 concludes (l45l l is convex regarding J\, therefore (|45T > is also convex regarding C?r, i.e. regarding 

{Gi}U " " m 

Lemma 9: The feasible set d46b for MIMO relay networks is convex regarding {G;}^. 

Proof: [36, §2.3] indicates that affine functions preserve convexity of sets, or allow us to construct convex sets from others. 
S C R" is a convex set. If / : R fe -)• R" is affine, i.e. /(x) = Ax + 6, where v4 G R nxfe and b G R n , the inverse image of S 
under /: f^ 1 (S) = {x\f(x) G S*}, is convex. 

Similar to Lemma [8] extract elements from {G/} as J-r, 5r in real domain. According to d48l ) and d42l . = Wr^r, 

where J-jg is composed of elements from {Uj }, {V^}, {Hj r, }, {HR l ] i }. [21] concludes that feasible set of J 7 ^ under constraint 
( f46b is convex, so that the feasible set of 5r as the inverse image of 6>r is also convex, i.e. the feasible set of {Gi}[ =1 is 
convex. 

■ 

Iterative RCRM algorithm for MIMO relay networks: Since lemma [8] and lemma [9] prove that (l47b is a typical convex 
optimization regarding {G;} for MIMO relay network, as well as regarding {Uj} and {Vi}. Then the optimization problem 
could be proceeded in an iterative manner, regarding {Vi}, {G;}, {Uj}, until convergence. At every step, the objective function 
(l45l l is non-increasing so that convergence is guaranteed. Compared with the algorithms relating to mean squared error (MSE) 
and sum rate as in J9], this algorithm directly optimizes DoF for high SNR for the MIMO relay network. 

2) MIMO Two-Way Relay Channel: Consider the if -pair two-way (bidirectional) relay MIMO interference channel as in 
11381 . The model is also as shown in Fig. [2] All Ik and Jk nodes are both sources and destinations, i.e. each Ik exchanges data 
with its partner node Jk via intermediate relaying nodes Ri, i?2, . . . , Rl- The whole procedure are three steps: uplink, forward, 
downlink. In uplink step, all end nodes Ik and Jk simultaneously transmit data towards relay nodes Ri, R2, . . . , Rl', in forward 
step, relay node Ri filters its received signal through a forwarding filter G;; in downlink step, relay nodes R\,Ri, . . . , Rl 
broadcast network-coded combined signal to all end nodes Ik and Jk- Denote H^^, H/^, Hj R t and H^j as channels 
from Ii to R h R t to I h R t to Jj and J 3 to R u respectively. R Rlh G C M * xM «, R hRl € C M e xM *, Ej jRl G C M - xM *, 
Hfl, jj G C Mx x Mu > . Although it does not affect our scheme, notice the downlink channel is reciprocal to the uplink channel, 
i.e. R RlIi = H IiRl , H Jji?! = R RlJj . 



Then the received signals Yj k and Y/ fc at Jk and Ik after self interference cancellation (also ignore noise term) are presented 

as: 

L K K 

1=1 i=i j=U&k 

K K 

= $ JkIk v Ik x Ik + J2 *^ i V/ i x /i + 2 , Vj X, 

i=l,i^k j — ^-^ji^k 

L L 

where $j kh = Y R, hRl GiR RlIi and $ JfcJ . = Y Hj^G/Hr, j. 
i=i 1=1 

L K K J 

Y/ fc =^H /fc H,Q( Hz.' / V/ X, • ^11,,., V, X., ) 

Z=l i=l,i^fe j'=l 

if if 

= $/ fc ./ fc V Jfc X Jfc + £ <i> IkU V Ii X Ii + Y: '"/ V / X -/ 

L L 

where ^ H/^jGjH^^ and ^H/^GjH^jj 

z=i z=i 

d49l shows the model is equivalent to 2K users transmitting in pairs, and there are (2K) 2 equivalent channels between 
transmitters and receivers which are {^j k i i }, {&.j k .j }, {®i k ii}, {®i k .J }■ The coding procedure is just the same as the case 
of MIMO relay networks. We could obtain the same RCRM problem for MIMO two-way relay networks as in d47i i. which is 
proved to be a convex optimization through the following similar lemmas. 

Lemma 10: For MIMO two-way relay networks as in d49l ), the corresponding rank minimization objective function d45l ) is 
convex regarding {Gi}f =l . 

Proof: Proof is similar to lemma [8] ■ 
Lemma 11: For MIMO two-way relay networks as in d49i >, the corresponding rank constrained feasible set d46*b is convex 
regarding {GJf =1 . 

Proof: Proof is similar to lemma [9] ■ 
3) MIMO Y Channel: Consider MIMO Y channel as in [39|. The model is also illustrated by Fig. [2] The network only has 
three user nodes which are set as I\, I2, h, and one relay node set as R\. Each user intends to convey independent messages 
for the other two users via the relay, while receiving two independent messages from the other two users. The transmission 
has two phases. First phase is the multiple-access channel (MAC) phase, when all users transmit to the relay; second phase is 
the broadcast (BC) phase, when the relay generates new transmitting signals and send to all users. Denote V^ 1 and X^ fcl as 
the precoding matrix and input data from 7^ to Ik- 

Then the received signals Yj k at Ik after self interference cancellation and Y Rl at R\ (both ignore noise terms) are: 



Y Rl =H Rl7l (V^X^+V^xfl) 
+ H flll3 (vK ll X^+V^X^]) 



(50) 

so that: 

Y 7l = (Hj lfil GiH iJl / 2 )V^ ll X^ 1 + (Hj^GiHr^V^ 1 X^ l] 

Y /2 = (H^GxH^M^X^ 1 + (H/^GxHr^JV^X^ 

Y /3 = (H/, Kl GiH JilJl )vg ,] X™ + (H^GxHr^V^X^ 

shows the model is equivalent to 6 pairs of users transmitting on 36 equivalent channels. We could obtain the same 
RCRM problem as in ( f4Tb for MIMO Y channel. It is proved to be convex through the following similar lemmas. 

Lemma 12: For the MIMO Y channel as in d50i >. the corresponding rank minimization objective function d45l l is convex 
regarding {Gi}f =1 . 

Proof: Proof is similar to lemma [8] ■ 



Lemma 13: For the MIMO Y channel as in d50l l, the corresponding rank constrained feasible set d46l ) is convex regarding 

Proof: Proof is similar to lemma [9] ■ 
Compared with [39], the analytical signal space alignment and network-coding-aware interference nulling beamforming 
scheme is complicated to be implemented and hard to adapt to changes of structure. While this RCRM method is a direct 
solution to form efficient and robust coding. 

4) MIMO Multi-hop Relay Channel: Consider MIMO interference channels with multi-hop relays as in Q, lfl8ll . The model 
in Fig. |2] needs to be modified. There are M layers of parallel relaying nodes between sources and destinations in the network. 
Denote the m-th layer's nodes as i?i, m , R2.n1, ■ ■ ■ , Rh,m, and the corresponding filter matrices are Gi, m , G2, m , . . . , G/, )TO . The 
channels are denoted as Hj.jj, M from M-th layer of relays to destination Jj, H^, m R l m l from (m — l)-th layer of relays to 
m-th layer of relays, 1 j i from source 7j to 1-st layer of relays. 

For simplicity, only assume there are two layers of relays. G„ 1 is the filter at q-th relay in first layer and G p 2 is the filter 
at p-th relay in second layer. Hr 2 r lt and Hj r 2 are channels from source Ii to relay Rq \, from relay R q> x to 

relay R p ,2, and from relay R p .2 to destination Jj. Then the received signal Yj j at Jj is: 

L L K 

Y ./j = / jHjjfl p i3 Gp |2 /j^-Rp ,2-Rg .lfig.i / ^Hfi g .aiVjiXj^ (51) 

p— 1 g— 1 i— 1 

The equivalent channel from Ii to Jj is: 

= Z)p=i Hj 3 -fl Pi2 G Pj2 X) g=1 H flp 2 fl 9 l G 9 ^H^^/, 

Iterative algorithms layer by layer: At each step, choose one layer of relays as the coding variables, then we could obtain 
the same RCRM problem as in (|47| i. The whole algorithm is convergent, since the objective function maintains the same all 
the time. It is proved to be a convex optimization problem through following similar lemmas. 

Lemma 14: For the MIMO interference channels with multi-hop relays as in (IBTT l. the corresponding rank minimization 
objective function (PES! is convex regarding {Gj,m}[f^wi IV 

Proof: Proof is similar to lemma [8] ■ 
Lemma 15: For the MIMO interference channels with multi-hop relays as in ( TSTI ). the corresponding rank constrained feasible 
set (06]) is convex regarding {G;, m }[^f:L (M) . 

Proof: Proof is similar to lemma |5j ■ 



In summary, MIMO interference channels with multi-hop relays are generally hard to have either analytical or numerical 
solutions 0, lfl8ll . While this novel RCRM based algorithm provides a robust and universal solution. 

5) Numerical Results for Diverse MIMO Relay Networks: All the numerical results are shown in Fig. [5] The above four 
scenarios of MIMO relay channel, MIMO two-way relay channel, MIMO Y channel and MIMO multi-hop channel are denoted 
as case 1, case 2, case 3, and case 4, respectively. Configure the four cases as shown in Table Hill In case 3, '1x2' means 2 
separate streams at one user in the MIMO Y channel. In case 4, '3x2' means 2 layers of relays (3 relay in each layer) in the 
multi-hop relay network. 
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TABLE III 

Configurations of 4 cases of MIMO relay networks 



Simulation has been done with the software tool CVX |[37l . In Fig. [5] 'Ref means random coding as a reference for the 4 
kinds of MIMO relay networks, which obviously obtain DoF. While 'Opt' means results with the IA coding design, which 
shows that DoF is well obtained. So that the novel application of RCRM method in relay networks successfully implements 
interference alignment in typical scenarios. Case 2 has the largest DoF because it has bidirectional transmissions. In case 3, 
MIMO Y channel transmits 2 datastreams for each user, however there is only 1 DoF for each virtual port. In case 4, multi-hop 
relays relaxe the requirements for individual relays in each layer to maitain high DoF. In summary, this novel relay coding 
design successfully implements interference alignment in above relay-involved scenarios as a effective universal solution. 




Fig. 5. Rates of 4 cases of MIMO relay networks (4 cases as configured in Table III before and after relay coding optimization) 



VII. Conclusion 

This work considers two strategies in which a finite number of relays can be employed to achieve IA in a i-C -pair interference 
network. In the coding-at-the-edge strategy, in the single antenna constant channel, relays can be used with time extension 
techniques to produce equivalent channels where standard IA schemes are applicable. However, at least two relays are proved 
to be required to produce channel randomness, which is a necessity for all kinds of IA schemes to be achievable. Novel 
precoding schemes with double-layered symbol extensions could also be used to investigate the issue of constant channels. 
In the coding-in-the-middle strategy, as a universal numerical approach for diverse general MIMO networks with relays, relay 
coding can be designed as an optimization problem to minimize the rank of the interference subspace. Efficient and robust 
algorithms are proposed. 
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